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Abstract 

Tidal torques play a key role in rotational dynamics of celestial bodies. They govern 
these bodies' tidal despinning, and also participate in the subtle process of entrapment 
of these bodies into spin-orbit resonances. This makes tidal torques directly relevant to 
the studies of habitability of planets and their moons. 

Our work begins with an explanation of how friction and lagging should be built 
into the theory of bodily tides. Although much of this material can be found in various 
publications, a short but self-consistent summary on the topic has been lacking in the 
hitherto literature, and we are filling the gap. 

After these preparations, we address a popular concise formula for the tidal torque, 
which is often used in the literature, for planets or stars. We explain why the derivation 
of this expression, offered in the paper by Goldreich (1966; AJ 71, 1-7) and in the 
books by Kaula (1968, eqn. 4.5.29), and Murray & Dermott (1999, eqn. 4.159), implicitly 
sets the time lag to be frequency independent. Accordingly, the ensuing expression for 
the torque can be applied only to bodies having a very special (and very hypothetical) 
rheology which makes the time lag frequency independent, i.e, the same for all Fourier 
modes in the spectrum of tide. This expression for the torque should not be used for 
bodies of other rheologies. Specifically the expression cannot be combined with an extra 
assertion of the geometric lag being constant, because at finite eccentricities the said 
assumption is incompatible with the constant-time-lag condition. 



1 



1 Context, Motivation, and Plan 

The mills of God grind slowly. 



Usually extremely weak, tidal interactions act upon celestial bodies for extended spans 
of time (up to billions of years). Over asons, tides shape celestial bodies' spin modes, and 
govern the exchange of angular momentum. The numerous and diverse manifestations of 
bodily tides range from the expected fall of Phobos to synchronous locking of the Moon, to 
Mercury's capture in the 3:2 spin-orbit resonance, to bloated hot-jupiter exoplanets in tight 
orbits around their host stars, to dissipational coalescence of short-period binary stars. This 
makes studies of tides essential to our understanding of the dynamical properties and evolution 
of stellar systems. 

The slow work of tides is responsible for circularisation, obliquity evolution, and synchro- 
nisation of planets and moons. The wide scope of these dynamical phenomena is not always 
matched by sophistication or versatility of the tidal models employed to describe them. 

1.1 Requirements to a consistent theory, 
and the history of simpler approaches 

A bona fide theory of bodily tides implies (1) decomposition of the tide into Fourier harmonic 
modes and (2) endowment of each separate Fourier mode with a phase delay and a magnitude 
decrease of its own. The first part of this programme, Fourier decomposition, was accomplished 
in full by Kaula (1964), though a partial sum of the Fourier series was developed yet by Darwin 
(1879). The second part of this programme, the quest for an adequate frequency dependence 
of the phase lags and dynamical Love numbers, is now in progress. While earlier attempts 
seldom went beyond the Maxwell model, more realistic rheologies are now coming into use. A 
rheology combining the Andrade model at higher frequencies with the Maxwell model at the 
lowest frequencies was investigated by Efroimsky (2012a,b)Q The necessity for such a combined 
model originates from that fact that different physical mechanisms of friction dominate tidal 
dissipation at different frequencies. Several other rheological laws were probed by Henning et 
al. (2009) and Nimmo et al. (2012). 

Some authors tried to sidestep a Fourier decomposition by building simpler toy models 
which would preserve qualitative features of the consistent tidal theory and, ideally, would 
yield some reasonable quantitative estimates. Two radically simplistic ad hoc tools, known as 
the constant geometric lag model (MacDonald 1964, Goldreich 1966, Murray & Dermott 1999) 

1 At lower frequencies, dissipation is predominantly viscous, and the mantle is well approximated with the 
Maxwell body. Its behaviour can be represented with a viscous damper and an elastic spring connected in 
series. Experiencing the same force, these elements have their elongations summed up. This illustrates how 
the total strain is comprised by a sum of viscous and clastic contributions generated by the same stress. 

At higher frequencies, the strain acquires the third component, one intended to describe inelasticity Inelas- 
ticity is produced by defect unpinning, a process effective at frequencies higher than a certain threshold (about 
1 yr _1 , in the case of Earth's mantle, - see Karato & Spetzler 1990). Hence, at frequencies above the threshold, 
dissipation is predominantly inelastic, and the mantle behaves as the Andrade body. 

A combined rheological model written down in Efroimsky (2012a, b) embraces both frequency bands. 

In our previous paper (Makarov et al. 2012), we made a slip of the pen by calling the Andrade creep anelastic. 
In fact, it should be called inelastic, which means: irrecoverable. (The term anelastic is applied to recoverable 
deformations, like the Maxwell behaviour.) 
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and the constant time lag model (Singer 1966; Mignard 1979, 1980, 1981; Heller et al. 2011; 
Hut 1981) are often resorted to, and are applied to rocky moons and planets and gas giants 
and stars alike. Historically, both these approaches were introduced for the ease of analytical 
treatment rather than on sound physical principles. 

1.2 Plan 

The first of the afore-said approaches, the constant geometric lag model , will be addressed in 
this paper. Our goal is to demonstrate that the model should be discarded, both for physical 
and mathematical reasons. On the one hand, the model, is not well grounded in the physical 
reality, because it assumes a constant tidal response independent of the rotation frequency 
everywhere except at the 1:1 resonance where it singularly changes sign. On the other hand, 
the model is genuinely contradictive in its mathematics. The source of the inherent conflict 
is the popular formula for the tidal torque (and its analogue for the tidal potential) through 
which the model is implemented. It turns out that these formulae tacitly imply Constance 
(frequency-independence) of the time lag, a circumstance prohibiting the additional imposition 
of the constant-geometric-lag ( = frequency-independent quality factor) condition. 

In the subsequent publication (Makarov & Efroimsky 2012), we shall explore the constant 
time lag model. This is an approach implying that all the tidal strain modes should experience 
the same temporal delay relative to the appropriate modes comprising the tidal stress. While 
the method works in the purely viscous limit (which, hypothetically, may be the case of stars 
and gaseous planets - see Hut 1981 and Eggleton 1998), there is no justification for using it for 
terrestrial-type bodies such as the Moon, Phobos, or any exoplanet with a rocky or partially 
molten mantle. In the light of the current rheological knowledge, the tidal response is very 
different, and its frequency dependence experiences especially steep variations in the vicinity of 
spin-orbit resonances. In Makarov & Efroimsky (2012), we shall demonstrate that illegitimate 
application of the constant time lag model to telluric objects leads to nonexistent phenomena 
like pseudorotation - not to mention that it squeezes the tidal-evolution timescales (Efroimsky 
& Lainey 2007) and alters the probabilities of capture into resonances (Makarov, Berghea, & 
Efroimsky 2012). 

2 The constant-torque model 

As we mentioned above, some authors tried to circumvent a consistent but laborious treatment 
by building simpler toy models. One such attempt was undertaken by MacDonald (1964, eqn. 
20) who assumed that the dynamical tide mimics a static tide, except for being displaced by a 
geometric lag angle. The same idea underlay the study, by Goldreich (1966), of a satellite on 
approach to the 1:1 spin-orbit resonance. Thus, for mathematical convenience, both authors 
set the geometric lag angle to be a frequency-independent constant. 

The approach is known as the constant angular lag model. This name, however, is inexact 
in the sense that, within the vicinity of the 1:1 spin-orbit resonance, the lags and torque 
change their signs twice over a period: when the bulge falls behind or advances (relative to 

2 Aside from its mathematical simplicity, the constant time lag method was motivated by its analogy with 
the viscously damped harmonic oscillator. This analogy, however, appeared in the literature a posteriori, 
Alexander (1973) being the earliest known to us work where this analogy is spelled out. 
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the direction towards the perturber), the sign is positive or negative, correspondingly. So in 
this discussion the term constant should be understood as frequency-independent: both the 
instantaneous phase lag and the instantaneous torque T are independent of the tidal frequency 
X- Consequently, the orbit-average (secular) tidal torque (T) is also frequency- independent. 

Sometimes this approach shows up in the literature under the name of MacDonald torque 
(e.g., Touma & Wisdom 1994, section 2.7.1). 

Unfortunately, the approach turns out to be inconsistent and should be discarded. Phys- 
ically, the constant angular lag model looks suspicious from the beginning, because in the 
vicinity of the 1:1 spin-orbit resonance it permits for abrupt switches of the torque, i.e., for 
situations where the torque changes its sign, retaining the absolute value. Although the abrupt 
switch can be substituted, by hand, with a continuous transition, this ad hoc alteration still 
would not save the method, because it would not heal a more fundamental defect. Mathe- 
matically, the derivation of the formula for the tidal torque within the said model contains a 
subtle and often unappreciated detail: this derivation implicitly sets the time lag At to be 
constant, as we shall demonstrate below. However, it can be shown that the assertion of the 
time lag being frequency independent is incompatible with the assertion of the geometric lag 
angle being frequency independent. In this way, the discussed simplified approach is inherently 
contradictive. 

Another defect of this approach is that it employs such entities as the instantaneous phase 
lag and the instantaneous quality factor, the latter being introduced as the inverse sine of the 
former. The so-defined instantaneous quality factor is not guaranteed to be related to the 
energy damping rate in the same manner as the regular (appropriate to a fixed frequency) 
quality factor is related to the dissipation rate at that frequency (Williams & Efroimsky 2012). 
Were the quality factor frequency-insensitive, this would not be a problem. However, the latter 
option is excluded within the discussed modeljf] as being incompatible with the constant-At 
assumption tacitly present. 

3 Goldreich (1966), Kaula (1968), Murray & Dermott 
(1999) 

Numerous authors offer the following expression for the polar component of the torque where- 
with the tide-raising perturber acts on the tidally-perturbed body: 

% = ^ G M? ^ k 2 sm2e g , (1) 

R being the radius of the disturbed body, Mi standing for the mass of the tide-raising per- 
turber, r denoting the instantaneous distance between the bodies, e g standing for the angular 

3 Rejection of the constant angular lag model does not, by itself, prohibit setting the quality factor constant, 
at least over some limited interval of frequencies. While realistic mantles never behave like this, such a rheology, 
in principle, is not impossible. However, employment of this rheology will not be available within a simplified 
model. Instead, one will have to attribute the same value to ki/Qi at all tidal modes, and then will have to 
insert this value of ki/Qi into all terms of the Fourier expansion of the tidal torque (the Darwin-Kaula series). 
Each such term will change its sign when the corresponding resonance gets transcended. 

Up to the late 60s of the past century, there was a consensus in the geophysical community that the seismic 
Q of rocks should be "substantially independent of frequency" (Knopoff 1964). This viewpoint was later 
disproved by a large volume of experimental data obtained both in the laboratory and in the field (see, e.g., 
Karato 2007 and references therein). 
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lag, and the obliquity being set nil. The subscript z serves the purpose of emphasising that 
the above formula furnishes the torque component orthogonal to the equator of the tidally 
perturbed body. 

Goldreich (1966) denotes the angular lag with \I> , Kaula (1968, eqn. 4.5.29) calls it 5, 
while Murray & Dermott (1999, eqn. 4.159) use the letter e. We shall follow the latter 
notation, though equipping it with the subscript g which means: geometric . 

Below we shall provide a detailed derivation of this formula, and shall see that the angle 
standing in it is indeed the instantaneous geometric tidal lag angle 



i.e., the instantaneous angular separation between the direction towards the bulge and that 
towards the perturber. Here At is the time lag, v is the true anomaly of the perturber, 9 is 

the sidereal angle of the perturbed body, and 9 is its spin rate. It is important to distinguish 
the instantaneous geometric lag e g from the instantaneous phase lag 



sometimes used in the literature (Efroimsky & Williams 2009, Williams & Efroimsky 2012). 

In Section 5 of Murray & Dermott (1999, eqns. 5.2 - 5.3), the authors rewrite the expression 
for the torque, employing (in fact, implying) the above expression of the lag through the true 
anomaly and spin rate: 

sin2e 3 = sin |2e 3 | Sgn (e g ) = - sin \2(9 - z>) At] Sgn(# - v) = - Q; 1 Sgn ( rj - up ) . (4) 

The new angles showing up in this formula are rj = 9 — Ai and ip = v — M. , with M. 
denoting the mean anomaly. These angles are depicted in Figure 5.1b in Ibid. Clearly, the 
time derivatives rj = 9 — n and ip = v — n are the spin rate and true anomaly rate in a frame 
which is centered on the tidally perturbed body @ and is rotating with the mean motion n . 
Interpreting the quantity Q s = 1/ sin |2(6 I — v) At\ as the instantaneous quality factor, the 
authors assume that they can set it frequency-independent (thus making the torque frequency- 
independent). 

This approach contains two flaws. First, as explained in Williams & Efroimsky (2012), it is 
not apparently evident whether the instantaneous quality factor introduced as the inverse sine 
of the instantaneous phase lag has the physical meaning usually instilled in a tidal dissipation 
factor at a certain sinusoidal mode. Second, and most important, is that in reality it is At 
which gets implicitly set as frequency independent in the derivation of ([T|). It then becomes 
impossible to assume that the geometric lag e g also is frequency-independent - the two as- 
sumptions are incompatible, as we shall see shortly. Consequently, setting the factor Q s to be 
frequency independent is no longer an option. This makes the entire constant geometric lag 
model or, to be exact, its implementation with ([[]), inherently contradictive. Specifically, it is 
illegitimate to assert that the tidal torque is proportional to Sgn(z> — 9) . 

4 In Murray & Dermott (1999), the role of a tidally perturbed body is played by the satellite, the planet 
acting as its tide-raising perturber. In a different setting, perturbed is the planet, the star or a satellite being 
the perturber. 



{i/ -6) At , 



(2) 



e ph = 2 (z> - 9) At = 2e g 



(3) 
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4 Mathematical Introduction. 



The tide-raising potential W created by the perturber always changes the shape and, as a 
result, the potential of the perturbed body. At the point R of the perturbed body's surface, 
the potential W(R, r *) created by the perturber can be expanded into a sum of terms 
Wi(R, r*) proportional to the Legendre polynomials P;(cos7) . Here 7 is the angle between 
the vectors r * and R pointing from the perturbed body's centre towards the perturber and 
the point on the perturbed body's surface, where the potential W is measured. 

Tidal distortion of the body's geometric form renders an addition U(r) to the body's 
potential in an exterior point r. This addition turns out to be comprised of terms Ui(r) each 
of which is proportional to the term Wi(R, r ) , with the surface point R located exactly 
below (i.e., having the same latitude and longitude as) the exterior point r* . 

The goal of this section is to provide a squeezed introduction into this formalism and to 
explain how it should be generalised from a static configuration setting to a dynamical setting. 

4.1 Static tides 

Let a body of radius R experience tides from a perturber of mass Ml placed at r* = 
(r*, <p\ A*) , with r* >R. 

— * 

At a point R = (R, 0, A) on the perturbed body's surface, the potential W due to the 
perturber is expanded over the Legendre polynomials P;(cos7) as0 

1 

W(R, r*) = W{H, r *) = - ^4 £ (4) P '( cos ^) 

1=2 T 1=2 \ 7 ' 

G M* ^ / R \ l ' (Z-m)> . . D ,. ,, D / . 

= -—^2^{ — ) 2^ ^ ! (2-^o m )Pm(sin0)P m (sm^) cosm(A-A*) 

1=2 ^ ^ m=0 ^ >' 

where 5ij is the Kronecker delta symbol, G is Newton's gravity constant, while 7 denotes the 
angular separation between the vectors r * and -R pointing from the centre of the perturbed 
body. The longitudes A, A* are reckoned from a fixed meridian on the perturbed body, the 
latitudes 0, 0* being reckoned from the equator. The integers / and m are called the degree 
and order, accordingly. The associated Legendre functions Pi m (x) are referred to as the 
associated Legendre polynomials when their argument is sine or cosine of some angle. 

The I th term Wi(R, r*) of the perturber's potential introduces a distortion into the 
perturbed body's shape, assumed to be linear. Then the ensuing I th amendment U\ to the 
perturbed body's potential will also be linear in Wi . Since Ui(r) falls off outside the body as 
r -(i+i) ^ overa n change in the exterior potential of the perturbed body will be: 

Uif) = = Z> ( ~ ) W t (R, r*) , (6) 

1=2 1=2 ^ ' 

5 Summation in formula (|5|) goes over I > 2 . The central term (I = ) is regarded as the principal, 
Newtonian, part of the potential generated by the outer body, and not as a part of the perturbation W caused 
by the finite size of the tidally perturbed body - indeed, the I = terms bears no dependence upon R . The 
reason why the I = terms is related to the fact that we are developing our formalism in the frame of the 
tidally perturbed body, not in an inertial frame. See, e.g., Efroimsky & Williams (2009, cqns. 5 - 11). 
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where k\ are the static Love numbers, R is the mean equatorial radius of the perturbed body, 
while R = (R, 0, A) and r = (r, 0, A) are a surface point and an exterior point above it, 
respectively, so that r > R . 

Combining ([6]) with §5$), we arrive at a useful formula for the amendment to the potential 
of the tidally disturbed body: 

U(r) = -GM*J2 k i i+i i+i J2 i. m | (2-M^m(sin0)PU^0*) cosm(A - A*) . (7) 
1=2 r r* m=0 y l + m )- 

This is how the tidally generated change in the perturbed body's potential is "felt" at a point 
r. The change is expressed as a function of the spherical coordinates r = (r, 0, A) of this 
point and the spherical coordinates r* = (r*, 0*, A*) of the tide-raising body. The formula 
may be employed when we have two exterior bodies: if one such body, a perturber of mass M x * 
located at r * produces tides on the perturbed body, then the other exterior body, located at 
r , will experience a potential perturbation ((?]) due to these tides. 

By changing variables from the spherical coordinates r* = (r*, 0*, A*) and r — (r, 0, A) 
to the Keplerian coordinates r* = (a*, e*, i*, Q*, u*, Ai*) and r = (a, e, i, Q, u, Ai ) , one 
obtains a formula equivalent to ((7j): 

™--i>(?f^(S)i:8^<' 

1=2 v 7 v 7 m=0 v ' 



-S 0m )^2Fi mp (i*) 22 G ipg( e *)z2 F imh(i) 22 Glh ^ e ) C0S {{ v *m Pq - m ®*) ~ ( v imhj ~ m9)) 

p=0 q=—oo h=0 j= — oo 

where 

vL pq = {l-2p)u>* + (l-2p + q)M* + mQ* , (9) 

v lmhj = (I -2h)u + (/ -2h + j)M + mVt , (10) 

q and j being arbitrary integers, p and h beng arbitrary nonnegative integers, Fi mp (i) 
being the inclination functions, while Gi pq (e) being the eccentricity polynomials coinciding 

(-l-l),{l-2p) 

with the Hansen coefficients X,, „ , , (e) . Also mind that 6* is the same as 6 , which is 

(l-2p+q) y I 

the sidereal angle of the tidally perturbed body. Following Kaula (1964), we equip 9 with an 
asterisk, when it shows up in expressions corresponding to the tide-raising body. In expression 
([8]), the terms — m6* and — m6 cancel one another, wherefore their presence may seem 
redundant. We better keep them, though, for they will no longer cancel when lagging comes 
into play. 

Decomposition ([8]) was pioneered by Kaula (1961, 1964). However, its partial sum, with 
\q\, \j\ < 2, was derived yet by Darwin (1879). In modern notations, Darwin's work is 

discussed by Ferraz-Mello, Rodriguez & Hussmann (2008). E 

For our further developments, it would be important to emphasise that formulae ([7j) and 

([H]) are equivalent to one another, because the latter is obtained from the former simply by a 

change of variables. 

6 Be mindful that the convention on the meaning of notations r and r* in Ibid, is opposite to ours. 
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4.2 Dynamical tides with no friction 

Derived for a static tide, formulae (J7J) and flS} extend trivially to an elastic dynamical setting 
where the tide adjusts instantaneously to the changing position r*(t) of the perturber. 

The key point is that, to get U(r) at the point r at time t , we insert into ((7j) or (JBJ) the 
perturber's position r*(t) taken at that same time t, and not at an earlier time. Formulae 
(J7J) and (jSJ) stay equivalent to one another, and remain unchanged, except that the sidereal 
angle and the angular coordinates acquire a simultaneous time dependence: 9 becomes 6{t) , 
<fi* becomes <fi*(t) , A becomes X(t) , u* becomes oo*{t) , £1 becomes Q(t) , while Ai* becomes 
M*(t). 

Thus, to obtain U (r(t) ) , we take the values of all variables at time t , leaving no place for 
any lagging. This is possible only for an absolutely elastic, i.e., frictionless perturbed body. 

4.3 Tidal modes and forcing frequencies 

Let us now write down the modes over which the tidal disturbance of the body gets ex- 
panded. We begin with expression ((SJ) for the perturbing potential at a fixed point R = 
(R, 0, A) on the surface of the perturbed body. Using the technique developed by Kaula 
(1961, 1964), we change the coordinates of the tide-raising body from r* = (r*, 0*, A*) to 
r* = (a*, e*, i*, Q*, u*, Ai* ) . However, the location on the body's surface, where the distur- 

— *■ 

bance is observed, is still parameterised with its spherical coordinates R = (R, (ft, A) : 

1=2 v 7 m=0 v ; 

l co , w — m even 

- 5 0m ) Pim(sm(j)) F imp(n G ipi( e *)\ r H npq -m(\ + e*)) . 

p=0 g=-oo ^ J l-m odd 

Since R is the radius of the tidally perturbed body, and since the latitude and the longitude 
A are reckoned from the equator and a fixed meridian, correspondingly, then (ITT]) is just another 
expression for the perturbing potential at the fixed point (R, (p, A) of the body's surface. In 
( ITT]) , the expression in round brackets can be reshaped as 

v L pq ~ m ( x + 0*) = ui mm (t - t ) - m A , (12) 

where 

u lmpq = (I- 2p) Co* + (/ - 2p + q) n* + m (Cl* - 9*) . (13) 

Here n* = Ai * is the mean motion of the perturber, while to is the time of perigee passage 
wherefrom the mean anomaly A4 * of the perturber is reckoned. 

We see from (IT2"]) that the quantities ui mpq given by (TLB"]) are the Fourier modes over which 
the tidal perturbation (fTT]) is expanded. While these modes can be positive or negative, the 
physical forcing frequencies, 

(14) 

7 The orbital parameters a , e , i , a* , e* , i* acquire no time dependence, insofar as the apsidal and nodal 
precession remain the only permitted variations of the orbits. 
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at which the stress oscillates, are positive-definite. 

Having developed formulae ([8]- [11]), Kaula (1961, 1964) never stipulated |f| that the Fourier 
modes of the tide are given by (fTBl . Possibly, he was not interested in the frequency dependence 
of the phase or time lags. 

In Section 6 of his book, Lambeck (1980) explained some aspects of Kaula's theory. While 
Lambeck's equation (6.1.13b) indicates that Lambeck could be aware of how the Fourier modes 
look, he too never wrote down the formula for the modes explicitly. Perhaps, like Kaula, 
Lambeck had no interest in the frequency dependence of lags - he just introduced a time lag 
At , which in his developments was implicitly regarded frequency independent. 

While in the review by Efroimsky & Williams (2009) and in Efroimsky (2012a,b) the ex- 
pression for ujimpq was written down explicitly, its derivation was omitted. 

Therefore, in the literature of which we are aware, the formula for the Fourier modes either 
was implied tacitly or was employed with no proof. This was our motivation to derive it here 
in such detail. 

To conclude, at the point R = (R, 0, A) of the perturbed body, the perturbing potential 
is expressed through the tidal Fourier modes as: 



OO / „ N 1 I 



GM* ^ f R\ ^ (l-m 



. a* J ^ (I + m)\ 

1=2 v 7 m=0 v ' 

i co , s i- m even 

- <W) Pi m {sm<f)) ^ F im P {i*) G 'w( e *) \ s j n f ( w Jmw(* - *o) - m\) 

p=0 q=-oo ^ J l-m odd 

In an idealised situation, when the extended body is frictionless and its response is instanta- 
neous, we can employ the static formula ([6]), as explained in subsection 14.21 Combining that 
formula with expression (115j) . we see that the additional tidal potential generated by a perfectly 

— * 

elastic body at the point r— (r, <p, A) right above R will now read as: 



i oo c n i-m even 

- S 0m ) P im (sin0) ^2 Fi mp (i*) ^2 G ipq( e *) k i \ sin \ (uim Pq (t - t ) - m\ 

p=0 q=-oo ^ J l-m Odd 

It is due to the lack of the friction that the Impq term of ( !T6|) is in phase with the Impq 
term of Below we shall see that inclusion of friction into the picture renders a phase shift 
between these terms. It is also in anticipation of the discussion of friction that we placed the 
Love numbers inside the ^2 mpq sum in expression (TIB"]) . Mode- independent in the perfectly 
elastic case, the Love numbers may acquire dependence upon the Fourier modes ui mpq , because 
friction may mitigate amplitudes of distortion differently at different frequencies. 



The linear combination standing on the right-hand side of our formula (|13p appeared in the denominator 
of formulae (29) and (50) in Kaula (1961). However, Kaula did not mention that this combination is a Fourier 
mode of the tide. 
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5 Friction and Lagging 



In this section, we shall trace, step by step, how internal friction gets included into the tidal 
theory. In subsection 14.21 we made an observation that, for an absolutely elastic (frictionless) 
body, treatment of dynamical tides mimics that of static tides, except that all coordinates 
acquire time-dependence. Our next step will be to incorporate friction, and therefore lagging, 
into the picture. 

As a first step, we shall address a simplistic method implying that the coordinates of the 
tide-raising body (as seen in a frame corotating with the perturbed body) get shifted back in 
time by some fixed time lag At . Although implementations of this method into formulae (J7J) 
and (jHJ) look very different, they render results which in fact are equivalent - simply because 
(J7J) and (jHJ) are equivalent, and because the same procedure (shift by At) is performed on the 
quantities with asterisks in both these formulae. 

However, the difference in the mathematical form of these, equivalent, results also prompts 
a more consistent way of taking care of friction. This, more advanced, method will be imple- 
mentable only in formula (jHJ) and not in (J7J). The method is the one used by Kaula (1964). 
Since the explanation of the method in Ibid, was extremely concise, we shall elucidate it here 
in mode detail. 



5.1 A naive way of bringing in lagging 

Naively, dissipation and the ensuing lagging can be included into the picture by assuming 
that the exterior body located at point r at time t is subject not to the tidal potential 
created simultaneously by the perturber residing at r*(t) , but to the potential generated by 
the perturber lagging in time on its orbit by some At . Speaking loosely, the no-asterisk 
exterior body located at r(t) "feels" the tide given by (J7J) or (jHJ), as if the tide were generated 
by the asterisk perturber located on its orbit not at r*(t) but at r*(t — At) . 

Mathematically, this implies that the tidal potential U(f) at t must be calculated via (J7J) 
or (jHJ), by insertion of the no-asterisk coordinates taken at t , and the coordinates with asterisk 
taken at t — At . The naive strategy also implies that the Love numbers ki keep their static 
values, though this detail is seldom spelled out. 

This approach implemented, our formulae (J7J) and (JHJ) will acquire the following form: 

21+1 



U(r(t)) = GMt'Ekt—-* ^£^4(2 

U r{t) l+1 r*{t-At) 1+1 ^ (' + "»)' 
- 8 0m ) P lm (sin <f)(t))P lm (sin <f)*(t- At)) cosm[X(t) - A*(t - At)} . (17) 



9 Here we imply: on its orbit as seen from the perturbed extended body. The caveat is needed, since we are 
considering the physical reaction of the extended body, and thus are interested in the location of the perturber 
relative to its surface, and not to an inertial frame. It is for this reason that in equation (|17[) we employ the 
latitudes and longitudes defined in a frame corotating with the perturbed extended body. In an inertial frame, 
a shifting of the perturber back by At should then be accompanied by a shift of the orientation of the extended 
body back by the same At ; and this is why we have 9*{t — At) in equation 1[T8|). 
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and 



EM? 



£=2 



GM* ( R 
a 



n \ x - (/ — m)! 

m=0 v * 



I 



6 



p=0 q=—oo 



h=0 



j=-oo 



COS 



( [v* lmpq {t - At) - me\t -At)}- [v lmhj {t) - m6(t) ] ) . (18) 

Just as their static precursors ([7]) and ([8]), our dynamical formulae (fT7|) and ({TBI remain 
equivalent to one another. They still reflect a mere switch from the spherical to the Keplerian 
coordinates, except that now f is taken at time t , while r * is taken at t — At . 

The longitude reckoned from a fixed meridian on the perturbed body is expressed through 
the true anomaly v , the periapse oo and the node Q as 

\ = -e + n + uj + v + o(i 2 ) . 

In neglect of the nodal and the apsidal precession, and for a small obliquity, this results in 

A « -6 + v + 0(z 2 ) . 
In (TIT]) , the argument of cosine may now be written, in a linear approximation over At , as 



m [A(i) - A*(t- At)] 



= m ( A - X 



m(9 - v) At + 0(i'- 



X(t) - X*(t) + A* At 

In (fT8|) . the argument of cosine may be shaped, in a linear approximation over At , as 
[v* mpq (t - At) - m 6*{t -At)]- \v lmhj {t) - m 6{t) ] = 



K P9 0) - m6*(t)] - [v lmhj (t)-m6(t)} 



(20) 



mpq 



where 

^Impq 



At 



(I - 2p) co* + (I - 2p + q) n* + m (fT - ^* 



At 



— ^Impq At (21) 

10 iJn route from (fT7|) to (JT8J) , one changes not only a coordinate system but also a frame of reference. Our 
longitude A being reckoned from a meridian, the switch goes from the corotating coordinates (r, <f>, A) to the 
Kepler coordinates (a, e, i, fl, u>, M) defined in a frame comoving but not corotating with the perturbed body. 

Technically, one first substitutes A with A — 9 , where A = A + 9 is the longitude in the comoving (not coro- 
tating) frame. Then one can resort to the standard formulae connecting the spherical and Kepler coordinates 
in the same frame. The formulae apply not to (r, <fi, A) but to (r, <fi, A) , see Kaula (1961). Thus, the current 
spin rate 9(t) pops up in (Tl8|) due to the transition from a corotating frame to a comoving one. 

All said relates equally to both the spherical and Kepler coordinates with asterisks. So the delayed value 
9*{t — At) , too, emerges in (TT5)) due to the frame switch. Recall that 9* is the same spin rate as 9 , except that 
it gets equipped with an asterisk, when it stands in expressions corresponding to the perturber. Also recall 
that, within the described approach, we model friction by simply shifting the perturber back in time by At. 

Leaving the coordinates (r, 0, A) untouched, and changing only (r*, <f>*, A*) to (a*, e*, i*, SI*, w*, M*) , 
one arrives at (fTT|) . Applying this machinery also to the variables with no asterisk, one ends up with (TPS]) . 
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is the phase lag corresponding to the mode ui mpq . 

From here we observe that the above-chosen method of taking the tidal friction into account 
fixes the phase lags in a very specific way: through shifting the perturber back on its orbit by 
a fixed time At , we set the phase lags ( 121 j) to be proportional to this At . 

Our formulae ( II 7p and ( II 8p can be written in another, equivalent form: 

tt,*> ^ GM * ^ ( R \ +1 ( R \ v-* {l-m)l tn , , 

1=2 V 7 V 7 m=0 v ' 

i oo (• -,1-m even 

P/ m (sin0) ^ Fi mp (i*) ^ Gi pq (e*) fcj <^ I (ui mpq (t - At - t ) - m\) 

p=0 q=-oo ^ ' l — m Odd 

This form is analogous to (fT6]) . except for the time lag At , the same for each Fourier mode. 

Expression (l2"2"j) is equivalent to expression ( 117]) . and is obtained from it by a switch from the 
spherical coordinates r * = (r*, </>*, A*) to the Kepler elements f* = ( a*, e*, i*, Q*, u*, M* ) , 
with the variables r = (r, 0, A) kept unchanged. From the equivalence of expressions (122]) and 
( TT7|) . we observe that, after the same time lag is added into all terms of ( )22|) . all the so-shifted 
terms add up to the equilibrium bulge geometrically displaced in such a way as if it were a 
static bulge generated by the perturber at a slightly different time. No other rheology can make 
this claim because, more generally, the lag in each term in the expansion would correspond to 
its own increment in t . This tells us that by setting At frequency-independent we impose a 
highly restrictive rheological rule, obedience to which cannot be expected of realistic mantles. 

A more profound problem of this approach lies in the fact that it is illegitimate to introduce 
lags, keeping at the same time the Love numbers unchanged. Mitigation of the magnitude and 
lagging of the phase are inseparably connected, though the link becomes apparent only within 
a consistent approach based on the Fourier expansion of the tide and on employment of one 
or another rheological law. That law will then define both lagging in phase and reduction in 
magnitude. 

5.2 A consistent way of bringing in lagging (Kaula 1964) 

The above expression (12T|) for phase lags contains in itself an obvious hint on how a general- 
type rheology should be built into the tidal theory - to that end, one simply has to endow 
each mode oui mpq with a time lag Ati(ui mpq ) of its own. Another adjustment is the mode 
dependence of the Love numbers: k% = ki(ui mpq ) . Expression ( TTg]) will now become 

1=2 v 7 v 7 m=0 K 1 

I oo / oo 

p=0 g=— oo h=0 j=—oo 

ki{uim P q) cos ( [v* mpq (t) - md*{t) ] - [vimhj(t) - mO{t) ] - ei(ui mpq ) ) } (23) 
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where 



ti{uim Pq ) = uj lmpq Ati(ui mpq ) = (I- 2p) uj* + (I - 2p + q) n* + m (fi* - 9* 



nipq J 



In (12"B1 - |2"3|) . we prefer to denote the phase and time lags not as ei mpq and Ati mpq , but as 
£i(wimpq) an d Atiiujimpq) ■ Indeed, their dependence on the indices mpq is solely due to the 
argument ui mpq . However, we cannot strip e or At of the subscript / , because the functional 
form of the frequency dependence of the phase or time lag is different for different / s. The 
tidal friction is not the same as the seismic friction, and the degree I affects the tidal damping 
rate0 Hence in formulae (1231 - [21]) we have e; and At; , and not just e or At. 

We would also write down the expression for the tidal potential in terms of the Keplerian 
elements of the perturber and the spherical coordinates of the point where this potential is 
observed: 



TTt* GM * ^ ( R \ +1 ( R \ (l-m)\ tn . , „ , . 



1=2 v ' v ' m=0 



i oo , w — m even 

Fim.p{i*) G>q(e*) ki(u)i mp q) < > (w/ mpg [t — At[(uJi mpq ) — t ] — m\ t 

p=0 q=-oo L k J / -m odd 

This form is analogous to ffTB 7 ]) and except for two details. First, each mode 0Ji mpq now 
has a time lag At; = Ati(u)i mpq ) of its own. Likewise, the dynamical Love number at each 
Fourier mode is a function of this mode: k\ = k[(ui mpq ) . 

The reason why we wrote down U (r, r*) in the above form is that it immediately furnishes 
an expression for the geometric lag angle of an arbitrary Impq bulge: 

filmpq — lmP9 At;(cj/ mpg ) . (26) 

m 

For example, the geometric lag angle of the principal, semidiurnal bulge is #2200 = U}2 2°° At = 

(n — 6) At , where the time lag is taken at the appropriate, semidiurnal mode: At = At2(cu22oo) • 
It is customary to introduce the convention that the phase and time lags and the Love 
numbers are functions not of the tidal mode ui mpq but of the positively defined frequency 
Ximpq = I ui mpq \ . Simplifying some calculations, this convention makes it necessary to intro- 
duce, by hand, sign factors into the terms of the Fourier expansions of the tidal force and 
torque (Efroimsky 2012b). 

In practical applications, most important is the special case when the exterior body located 
at r coincides with the tide-raising perturber located at r * . In this situation, the perturber 
is acting upon itself through the tide it creates on the perturbed body. Keeping the phase 
lags intact, and identifying r(t) with r *(t) , we may be tempted to mis-assume that the 
expression with no asterisk, standing in the argument of the cosine in ( 12 3 j) . compensates the 
expression with an asterisk, so the phase lag becomes all that is left. However, this would only 
furnish us the secular term of the tidal potential U (r) wherewith the perturber acts upon itself 



11 The difference between the tidal and seismic friction and, accordingly, the difference of dissipation at 
different / s is unimportant in small bodies, where only the rheology matters. However things change in large 
planets where self-gravitation becomes a crucial factor in tidal friction (Efroimsky 2012a). 
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through the medium of the tidally perturbed body. This term is proportional to cos ei mpq . As 
the p, q and h, j are independent pairs of integers, there also will be contributions with 
{Pi o} {h, j} ■ These are the oscillating components of the tidal potential U(r) . Although 
their time average is nil, they do contribute to the heat production, while the appropriate 
oscillating components of the tidal torque may influence free librations. The topic was first 
addressed in Efroimsky (2012b), and was revisited by Makarov et al. (2012) who explored 
whether the oscillating part of the torque can influence capture into spin-orbit resonances. 
It has turned out that, naturally, the oscillating part of the torque alters the outcome of a 
particular realisation of the capture scenario, but leaves the statistics unchanged. 

6 Calculation of the tidal torque 

We would begin with emphasising a key circumstance concerning the two forms of expansion 
of the tidal potential. In the absence of friction, these expansions, (J7J) and (jSJ), were equivalent. 
Their amended versions, ( 1T7|) and f[T8"j) . remain equivalent in the presence of friction, provided 
the latter obeys a special (often unphysical) restriction that the time lag At is the same at 
all frequencies. Beyond that threshold, the equivalence of the two expansions failed to extend. 
Since expression ( ITT]) does not contain the Fourier modes in it, it is plainly impossible to write 
(II 7p in a form that takes into account different time lags at different frequencies. Therefore, 
any calculation based on (j!7p will unavoidably imply the frequency-independence of At and 
will thus be incompatible with any other rheological law. 

6.1 Tidal torque, as derived from the concise expression (1171) 

Let us calculate the tidal torque, using (TT7]) . By employing this expression, we automatically 
set the rheology to be At = const . 

Consider an exterior body of mass Mi located at r, which is subject to the additional 
tidal potential U (r) of the tidally perturbed body. Then its energy in this potential will be 
Mi U {f) . When the position of the exterior body is rendered by the spherical coordinates, 
the polar component of the torque acting on it can be conveniently expressed as: T z = 
— Mi dU(r)/dX . The polar torque wherewith the exterior body acts back on the tidally 
perturbed body is the negative of T z : 



Here polar means: orthogonal to the perturbed body's equator. For small obliquities (and, 
therefore, small latitudes), insertion of (fT7| into (|2T|) yields 



r z {f) 



Mi 



i 



dU(r) 



(27) 



%(r(t)) 



- <W) Pi m (Q)Pim(Q) sin ( m (A - A*) - m {9 - v) At ) + 0(f) , (28) 




where we made use of (|T9|) . 
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Since the integer m is now entering the above expression as a multiplier, the term with 
m = becomes nil. As P 2 i(0) = 0, the term with m — 1 also vanishes. Hence, the term 
where / = m = 2 is leading. Neglecting the smaller terms, we thus obtain: 

T z {r{t)) « jj GMi M?k 2 — 3 sin (2 (A- A*) - 2 (9-0) At) . (29) 

2 r(t)V(t- At) V 7 

Consider the special case when the tide-raising perturber (with the asterisk) coincides with the 
other external body (with no asterisk). The perturber creates tides on the perturbed body, 
and then interacts with the tides it itself has created. Hence the perturber becomes subject to 
a tidal torque T exerted by it upon itself, through the medium of the tidal bulge it creates 
on the perturbed body. Evidently, a torque T = — T , of the same magnitude but opposite 
direction, will be acting upon the perturbed body. We arrive at this torque by setting A = A* 
and Mi = M* in the above expression: 

%(r) « ^GM^^sin^-^At) . (30) 

Naturally, for the model with a frequency-independent At , the quantity 

e g = (it- 6) At (31) 

is the geometric lag, i.e., the angular separation between the planetocentric directions towards 
the perturber and the bulge. Accordingly, within the said model, the quantity 

X = 2\v-6\ (32) 

acts as an instantaneous tidal frequency. 
The quantity 

e ph = 2 (v - 6) At = 2e g (33) 

is commonly interpreted as an instantaneous phase lag , so the torque in expression (130]) may 
be written down as 

T z (r) ^*-GM?k 2 \ sine p , = ^ G M, 2 k 2 ^ sin2e 5 , (34) 

which is exactly the expression ([T|) of our concern. In Section 2, we mentioned several popular 
papers and books, including Murray & Dermott (1999, eqn. 4.159), where this formula is 
employed. Now, that we have derived this formula accurately, we see that its validity hinges 
on the time lag being frequency independent. 



12 Be mindful that in ([50)1 we chose to make no distinction between r(i) and r(t — At) . Replacement of 
r(t — At) with r(t) gives birth to an absolute error of order 0(eQ~ 2 n/x) ■ It is however explained in Efroimsky 
& Williams (2009), that after averaging of (|30|) over one orbital period this error reduces to 0(e 2 Q~ 3 n 2 /x 2 ) ■ 

13 Mind a misprint in Eqn. (4.159) of Murray & Dermott (1999): in the denominator, a 6 must be changed 
to r 6 . The misprint emerged because in subsection 4.2 the distance was denoted with a . In formulae (5.2 - 
5.3) of Ibid, the misprint gets corrected. 
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While it is common (McDonald 1964, Goldreich 1966, Kaula 1968, Murray & Dermott 
1999) to treat 



Q = 1/| sin(2(z>-0) AH | = l/sm\e ph \ (35) 

as an instantaneous quality factor, the validity of this interpretation of ( 1351) remains question- 
able. For a nonzero eccentricity, the instantaneous tidal frequency ( 1321) is varying in time. So 
it is not readily apparent whether the instantaneous Q is connected to the damping rate in 
the manner the proper quality factor introduced at a certain frequency links to the damping 
rate at that frequency. McDonald (1964) and Goldreich (1966) tried to sidestep this difficulty 
by assuming that Q is a frequency independent constant. However, at finite eccentricities this 
assumption does not work, as it is incompatible with the constant At assumption (the latter 
assumption being a necessary prerequisite to using formulae I3T1 - 135| as we saw above). For 
more on this see Williams & Efroimsky (2012). 

6.2 Tidal torque, as derived from the Fourier expansion (1181). 
with all Fourier modes delayed by the same time lag At 

When starting out with expression (ITS]) , it is convenient to use the formula 

T.m = - Mi «m , (36) 

instead of ( 12~T1) . Technically, we should first differentiate U with respect to the sidereal angle 
9 , and then set 9 = 9* . We should also set the orbital elements with asterisk equal to their 
counterparts with no asterisk, it being understood that the tide-raising perturber is the same 
as the other exterior body which "feels" the tides on the perturbed body. The development will 
furnish us the polar component of the torque with which the perturber acts upon the tidally 
deformed body. For exploration of dynamics in a low-obliquity configuration, this component 
is sufficient. 

Insertion of the Fourier series f fTB"]) into equation f|36|) yields a Fourier series for the polar 
component of the torque, which is presented in Efroimsky (2012b). Here we shall not repeat 
this long formula, but shall only make an important comment on it. Insofar as At stays 
frequency independent (i.e., has the same value for all phase lags ( 12TT) entering the expansion 
for the torque), the resulting series for the torque stays fully equivalent to ( )28|) . with A and 
A* set equal to one another in the latter formula. 

This equivalence is ensured by the expression ( ITS]) for the potential U being equivalent 
to the expression (fT7|) whence formula (1281) originated, and by our agreement to keep At the 
same for all phase lags. 

6.3 Tidal torque, as derived from the Fourier expansion (1231). 
with each mode LUi mpq having a time lag of its own, Ati(ui mpq ) 

As soon as we abandon the assumption that the phase lags (12~T1) contain the same fixed At , i.e., 
as soon as we switch from the lags (12~T1) to those rendered by (12"4|) . we acquire an opportunity 
to describe a tidal torque acting on a perturbed body of an arbitrary rheology. Indeed, as the 
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time lags can have an arbitrary mode- dependence, this also relates to the phase lags. Above 
that, we now permit the Love numbers to be mode- dependent. 

To derive the tidal torque, we now combine formula (136]) not with expansion (11 8p but with 
the expansion ( 12 3 j) where the time lags are, generally, all different. The rheological emancipa- 
tion, though, comes at a cost: the Fourier decomposition for the torque, obtained through ( )2"3"|) . 
with mode- dependent Ati(u lmpq ) and ki(uJ lm ) , will no longer be equivalent to the concise 
and elegant formula ( 128|) . The customary and widely used leading-order approximation of ( 128]) . 
given by ( 130]) or by (134"]) . will not work for an arbitrary rheology. 

ki(ui, ) 

If, for example, we choose to set the factors jy> — mpq >. Sgn(o;, ) = ki sine; to be mode 

independent, then, to calculate the torque, we shall have to plug the same value of ki sine; 
into all terms of the Fourier series for the torque (eqn. 106 from Efroimsky 2012b). However, 
we shall not be able to employ the neat formula (I31E]) . 



7 Why the constant angular lag ' model is wrong 

As we explained above, an accurate derivation of the popular formula (|T|) for the polar com- 
ponent of the torque hinges on a tacit assumption that the time lag At is the same for all 
modes in the expansion of the tidal potential. Through formula (131 p . this, mode-independent 
time lag is related to the geometric lag angle e g present in formula (Q. 

Relation ( )3T|) tells us that, since At is the same for all the tidal modes involved, then 
the geometric lag angle e g cannot be treated as a fixed constant - even if we add a caveat 
permitting e g to switch signs when the value of v transcends 9 . With this caveat or not, 
keeping | e g | constant is impossible simply for the reason that (for a nonvanishing eccentricity) 
the quantity v oscillates in time 1^1 So a constant At is incompatible with a constant | e g \ . 

This is the reason why the so-called constant angular lag model based on ([T|) must be 
discarded wholesale as being inherently contradictive. 

k[ (u, ) 

While we still retain the right to set the factors jy-, — mpq >. Sgn(a; ;mp(J ) = k\ sine; mode 

independent, the value of the torque resulting from this assumption has to be calculated by 
insertion of these factors into the full Fourier expansion of the torque and not into ([I]). We 
may as well use (CQ), but only for a constant time lag, and not for a constant geometric lag 
angle. 



8 Conclusions 

We have reexamined the common formula ([I]) for the tidal torque, a formula which is equivalent 
to the expressions given in Sections 4 and 5 of Murray & Dermott (1999) and to the expressions 
offered in Goldreich (1966) and Kaula (1968). It has turned out that an accurate derivation 

14 While derivation of (fT]) absolutely requires At to be the same for all tidal modes, it does not require 
At to be fixed in time. Therefore, in theory, we can save the constant angular lag model by tuning the time 
dependence of At in such a special way that e g in f|3 If) stays constant in time. This however would require the 
dissipativc properties of the mantle to be fine-tuned, simultaneously at all frequencies, to ensure that, first, the 
time lags at all frequencies evolve but remain equal to one another and, second, that At [y — 8) stays constant 
in time. As the evolution rate of v — 9 is defined by the orbit, such fine tuning of rheology is unrealistic. 
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of this popular formula necessarily implies a specific rheology - the assertion that the time lag 
At is frequency independent. As can be easily seen from (1311) . this assertion is incompatible 
with the assertion of the geometric lag being frequency independent. Moreover, the quantity 
e g furnished by formula (1311) can be endowed with the meaning of a geometric lag only within 
the constant At rheological model (and only for a small obliquity i). 

To conclude, whenever the analysis of bodily tides is carried out using ([[]), the analysis can- 
not be combined with a constant geometric lag (or phase lag, or quality factor) assumption, 
nor with any other assumption different from the frequency independence of At . This circum- 
stance would not, by itself, prohibit one from considering a material for which the factors 
k\ {oJ lm ) 

tw — mpq \ Sgn(a;, ) = k\ sine; are insensitive to the frequency over some limited frequency 

band. Consistent employment of this model will then require insertion of the same value of 
ki sin ei into all terms of the expansion of the torque (each term thus changing its sign as the 
corresponding resonance is transcended). However, neither the quantity (z> — 9) nor its sign 
will come into play in this expression for the torque. So the outcome will be different from the 
mathematically incorrect "constant angular lag" model based on equation (Tj[|). 
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